Optical forces and torques in non-uniform beams of light 
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The spin angular momentum in an elliptically polarized beam of light plays several noteworthy 
roles in optical traps. It contributes to the linear momentum density in a non-uniform beam, and 
thus to the radiation pressure exerted on illuminated objects. It can be converted into orbital angular 
momentum, and thus can exert torques even on optically isotropic objects. Its curl, moreover, 
contributes to both forces and torques without spin-to-orbit conversion. We demonstrate these 
effects experimentally by tracking colloidal spheres diffusing in elliptically polarized optical tweezers. 
Clusters of spheres circulate determinisitically about the beam's axis. A single sphere, by contrast, 
undergoes stochastic Brownian vortex circulation that maps out the optical force field. 



Optical forces arising from the polarization and po- 
larization gradients in vector beams of light constitute 
a new frontier for optical micromanipulation. Linearly 
polarized light has been used to orient birefringent ob- 
jects in conventional optical tweezers [lHH and circular 
olarization has been used to make them rotate P, 0- 
More recently, optically isotropic objects also have 
been observed to circulate in circularly polarized optical 
traps . through a process described as spin-to-orbit 

conversion 

EMI- Here, we present a general formu- 
lation of the linear and angular momentum densities in 
vector beams of light that clarifies how the amplitude, 
phase and polarization profiles contribute to the forces 
and torques that such beams exert on illuminated ob- 
jects. This formulation reveals that the curl of the spin 
angular momentum can exert torques on illuminated ob- 
jects without contributing to the light's orbital angular 
momentum, and that this effect dominates spin-to-orbit 
conversion in circularly polarized optical tweezers. Pre- 
dicted properties of polarization-dependent optical forces 
are confirmed through observations of a previously un- 
reported mode of Brownian vortex circulation for an 
isotropic sphere in elliptically polarized optical tweezers. 

The vector potential describing a beam of light of an- 
gular frequency uj may be written as 



A(r,t)=u(r) e l ^~ lujt e{r), 



(1) 



where u{r) is the real-valued amplitude, <p(r) is the real- 
valued phase and e(r) is the complex- valued polarization 
vector at position r. This description is useful for prac- 
tical applications because u(r), tp(r) and e(r) may be 
specified independently, for example using holographic 
techniques Poynting's theorem then yields 

the time-averaged momentum density 



9( r ) = 



2pic 2 



'{A*(r,t)x [VxA(r,t)]}, (2) 



where \x is the permeability of the medium and c is the 
speed of light in the medium. The momentum density 
gives rise to the radiation pressure that the light exerts 
on illuminated objects and may be expressed in terms of 



the experimentally accessible parameters as 
ff(r) = j^/(r) V^-^/(r)e*Ve, + ivx S) (3) 
where I(r) = u 2 (r) is the intensity and where 



s(r) 
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I(r) ar(r) 



(4) 



is the spin angular momentum density in a beam of light 
with local helicity 



cr(r) = i e(r) x e*(r) 



(5) 



The projection of er(r) onto the propagation direction 
k(r) is related to the Stokes parameters of the beam [Uj 
by <r(r) ■ k(r) = Ss(r) /So(r). It achieves extremal val- 
ues of +1 and —1 for right- and left-circularly polarized 
light, respectively. 

The momentum density described by Eq. gives rise 
to the radiation pressure experienced by objects that ab- 
sorb or scatter light. Identifying g(r) with the radiation 
pressure on a particle is most appropriate in the Rayleigh 
limit, when the particle's size is no greater than the wave- 
length of light. In this limit, the three terms in g(r) may 
be interpreted as distinct mechanisms by which a beam 
of light exerts forces on illuminated objects. 

The first two terms in Eq. ([3]) constitute the famil- 
iar phase-gradient contribution to the radiation pressure 
[3]. In this context, the second term accounts for the 
independent phase profiles that may be imposed on the 
real and imaginary components of the polarization in an 
elliptically polarized beam. Phase gradients have been 
used to create three-dimensional optical force landscapes 
[lij ]. such as knotted force fields [l9| and true tractor 
beams [2fj| . They also account for the orbital angular 
momentum density 



£(r) 
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I(r) [r x {V<p - ie*Vtj)} , 



(6) 



carried by helical modes of light [21|, |22|. In this con- 



text, the polarization-dependent term in Eq. ([6]) vanishes 
identically in linearly polarized light, but manifests spin- 
to-orbit conversion in elliptically polarized beams. 
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FIG. 1: (Color online) (a) Streamlines of the momentum density g(r) in a right-circularly polarized optical tweezer, (b) 
Components of g(r) in the plane indicated in (a), shaded by the intensity I(r). (c) Measured trajectory of one particle in a 
seven-sphere cluster trapped near the focus of the beam. Discrete points show the last three seconds of motion, colored by time, 
(d) Circulation rate Q, as a function of the beam's Stokes parameters S3 /So- Inset: snapshot of the cluster indicating the sphere 
whose trajectory is plotted, (e) Three seconds of a 3.5-minute trajectory of a single polystyrene sphere diffusing in a circularly 
polarized optical tweezer, shaded by time, (f) Time-averaged probability flux j(r) computed from the full measured trajectory. 
Barbs are colored by the relative probability density p(r) computed from the same trajectory. Brownian vortex circulation is 
apparent in the vorticity of j(r). (g) Dependence of the Brownian vortex circulation rate on S3/S0. Inset: snapshot of the 
trapped sphere. The color bar indicates relative intensity 7(r) for (b), time for (c) and (e), and relative probability p(r) for (f). 



The third term in Eq. describes how variations in 
spin angular momentum contribute to the linear momen- 
tum density in non-uniform beams of light. This spin- 
curl term encompasses forces due to spatially- varying el- 
liptical polarization and also those due to intensity varia- 
tions in elliptically polarized beams. Streamlines of V x s 
naturally loop around extrema in the beam's intensity. 
Spin-curl forces thus tend to make illuminated objects 
circulate in the plane transverse to the direction of prop- 
agation. Observations of colloidal spheres circulating in 
beams of li ght with spatially- varying elliptical polariza- 
tion consequently have been interpreted as ev- 
idence that the curl of the polarization contributes to the 
light's orbital angular momentum. Equation ©, how- 
ever, makes clear that the spin-curl contribution to g(r) 
does not contribute in any way to £(r). For the same 
reason, observations of optically-induced circulation in 
uniformly circularly-polarized optical traps [H [t| need 
not imply spin-to-orbit conversion. 

To illustrate these point, we consider the forces ex- 
erted on an optically isotropic colloidal sphere by ellip- 
tically polarized optical tweezers. We model the trap as 
an Gaussian beam of wavenumber k brought to a focus 
with convergence angle a by a lens of focal length / and 
numerical aperture N A = n m sin a in a medium of re- 



fractive index n m . The beam's initial polarization is 



1 



e(r) = -j, (x + 



y 



(7) 



with a corresponding incident helicity <jq — sin 5 along z . 
The focused beam's vector potential may be expressed in 
cylindrical coordinates r = (p, <t>, z) with the Richards- 
Wolf integral formulation 0, 0, [23 , 



A{r ) = -i [A (r ) + A 2 (r )] (cos <p + e iS sin 4 
- i [A Q (r) - Ai (r)] (e lS cos <p - sin cj>) 
-2Ax(r) (cos(f> + e iS sin ( 
as a Fourier-Bessel expansion 

k f u o f° _ m j „«fe< 



A n (r) 



o. n {0) J n (kp sin 9) e 1 



2iu> jo 

with expansion coefficients (HJ 

oo(0) = (1 + cos 6) sin 0V cos 9 
a 1 (9) = sin 2 9 Vcos 9 



z, (8) 



d9, (9) 



012(6) = (1 — cos 9) sin 9 vcos9. 



(10) 

(11) 
(12) 



Streamlines of g(r) in a right-circularly polarized optical 
tweezer (ctq = +1, NA = 1.4) are shown spiraling around 
the optical axis in Fig. [lja). 
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A slice through the beam in the transverse plane in- 
dicated in Fig. HJa) reveals the azimuthal component to 
the transverse momentum density g±(r) = g(r) ■</> that 
is plotted in Fig. ([Ud). The transverse momentum den- 
sity may be resolved into two contributions 



9±{r) = go{r)+gs(r) 



(13) 



arising from the spin-to-orbit and spin-curl contributions 
to g(r), respectively: 

go(r) = ^-\\A 1 (r)\ 2 + \A 2 (r)\ 2 ]a and (14) 
/iC z p L J 

9s(r) = [d^{At(r) (A (r) - A 2 (r))} - 



(\A (r)\ 2 -\A 2 (r)f 



cr . 



(15) 



Both are proportional to the helicity of the incident 
beam, Co- They do not, however, contribute equally to 
the transverse component of the radiation pressure. At 
the focus of the circularly-polarized optical tweezer, for 
example, 79% of the transverse momentum density is due 
to the spin-curl term gs(f) and only 21% from spin-to- 
orbit conversion. More generally, both Ai(r) and A 2 (r) 
vanish in the paraxial approximation; there is no spin-to- 
orbit conversion in weakly focused beams. The spin-curl 
contribution, by contrast, persists in the paraxial limit. 

We probe the properties of spin-dependent optical 
forces by measuring their influence on the motion of 
micrometer-scale colloidal spheres. Our system con- 
sists of 1.0 |J.m diameter polystyrene (PS) spheres (Poly- 
sciences, Lot # 586632) dispersed in water and trapped 
in optical tweezers whose helicity Co is controlled with 
a quarter- wave plate. The isotropic dielectric spheres 
absorb very little light directly. By scattering light, how- 
ever, they experience radiation pressure proportional to 
the local momentum density. Our optical tweezer is pow- 
ered by up to 4 Wof laser light at a vacuum wavelength 
of A = 532 nm (Coherent Verdi 5W). The elliptically 
polarized beam is relayed with a dichroic mirror to the 
input pupil of an objective lens (Nikon Plan Apo, 100 x, 
NA 1.4), which focuses the light into a trap. We account 
for the mirror's influence on the polarization by mea- 
suring the beam's Stokes parameters in the input plane 
of the objective lens. The sample is imaged using the 
same lens in conventional bright-field illumination, which 
passes through the dichroic mirror to a video camera 
(NEC TI-324AII). Digitally recorded video is analyzed 
with standard methods of digital video microscopy [26| 
to measure the trajectory rj = r(jr) of a probe particle 
with 10 nmresolution at r = 33 ms intervals. 

The trajectory plotted in Fig. QJc) was obtained for 
one of seven spheres trapped against a glass surface by 
a right-circularly-polarized optical tweezer (cr = +0.8) 
powered by 1.5 W. The optically-assembled cluster, 
shown inset into Fig. [IJd), spans the region of the beam 



indicated in Fig. Hfb), and thus rotates about the beam 
axis at a rate of roughly il — 0.4 Hz. The data in 
Fig. [ljd) confirm the prediction of Eqs. (fT4|) and (fT5|) 
that the rotation rate varies linearly with the degree of 
circular polarization. 

The colloidal cluster circulates deterministically in the 
elliptically polarized optical tweezer because it continu- 
ously scatters light in regions where g±(r) is substantial. 
A single sphere diffusing in an elliptically polarized op- 
tical tweezer, by contrast, explores the entire force land- 
scape presented by the light. This includes regions near 
the optical axis where g±(v) is predicted to vanish. Fig- 
urelTfe) shows the measured trajectory of one such sphere 
in a right-circularly-polarized trap (uo = +0.8) powered 
by 0.05 W. Optically-induced circulation is not immedi- 
ately obvious in the noisy trajectory, which is shaded to 
indicate the passage of time. It becomes evident when 
the trajectory rj is compiled into a time-averaged esti- 
mate [27] for the steady-state probability current 



j(r) = 



1 



N-l 

\- r j+ i - rj 



N — 1 ^— ' 



r — 



(16) 

which is plotted in Fig. QTf). Here N = 7,000 is the 
number of discrete samples, and 8 a {r) is the kernel of 
an adaptive density estimator (27j whose width a varies 
with the sampling density. The symbols in Fig. [IJf ) are 
shaded by the estimated probability density 



P(r) 



N ^ 



N 



(17) 



for finding the particle near r. Together, j(r) and p(r) 
confirm the prediction of Eqs. (fT4"|) and (fT5"j) that circu- 
lation vanishes on the optical axis where the particle's 
probability density is greatest. 

Taking care to measure r from the center of circulation, 
the mean circulation rate may be estimated as 



Q, = / p(r) [r x j(r) 



■ zd 2 r. 



(18) 



Equation (|18[) improves upon the graphical method for 
estimating introduced in Ref. 1 28] by making optimal 
use of discretely sampled data [27J. Because the single 
particle spends most of its time in a curl-free region of 
the optical force field, its circulation rate is substantially 
smaller than in the deterministic case. Even so, the data 
in Fig. [IJg) again are consistent with the prediction that 
n scales linearly with <7o. 

The single particle's stochastic motion differs qualita- 
tively from the cluster's deterministic circulation. Were 
it not for random thermal forces, the isolated sphere 
would remain at mechanical equilibrium on the optical 
axis. Thermal forces enable it to explore the optical force 
landscape, where it is advected by the spin-dependent 
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contribution to the radiation pressure. This system, 
therefore constitutes an example of a Brownian vortex 
(29l l30j |. a stochastic machine that uses noise to trans- 
duce work out of a static non-conservative force field. 
Unlike previous experimental demonstrations of Brow- 



nian vortexes |28|, |29[ the conservative restoring force in 
this system is transverse to the non-conservative con- 
tributions. Consequently, the particle's radial excur- 
sions are described by the Boltzmann distribution [3(| 
p(r) = exp (— pU(rj) where /3 _1 = fcgT is the thermal 
energy scale at absolute temperature T and U(r) is the 
particle's potential energy in the trap. The solenoidal 
part of the radiation pressure [30j ] 
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f(r) 



87T/iC ; 



V x 



V'xJ(r') [VV-ieJV'ej] 



dV 



+ |vxa(r), (19) 



then advects p(t") into the probability current j(r) — 
H P p(r) f(r), where £ is the particle's scattering cross- 
section and \i v is its mobility. The first term in Eq. (fT9|) 
may be neglected in a beam such as a optical tweezer 
that carries little or no orbital angular momentum. The 
data in Fig. [ljf) thus map out the spin-curl force in the 
transverse plane. Moreover, because the circulation di- 
rection is determined unambiguously by the curl of f(r), 
this system is a practical realization of a so-called trivial 
Brownian vortex, which has been proposed [30j j but not 
previously demonstrated. 

Formulating the optical momentum density in terms 
of experimentally accessible parameters clarifies the na- 
ture and origin of the forces that can be applied to mi- 
croscopic objects using the radiation pressure in beams 
of light. This formulation confirms previous reports of 
forces arising from phase gradients [l6| and demonstrates 
that phase-gradient forces act independently of the state 
of polarization. The spin-curl mechanism unifies forces 
arising from the curl of the polarization and forces due 
to intensity gradients in elliptically polarized beams. Be- 
cause they induce circulatory motion, spin-curl forces are 
easily misinterpreted as evidence for spin-to-orbit con- 
version. The spin-curl density, however, does not con- 
tribute to the orbital angular momentum of the light. 
Spin-to-orbit conversion, by contrast, removes spin an- 
gular momentum from a beam o f lig ht and transmutes it 
into orbital angular momentum . The present formu- 
lation clarifies this mechanism, and reveals that spin-to- 
orbit conversion has played a secondary role in previous 
reports of optically-induced circulation. Using Eq. ([3]) 
as a guide, all three mechanisms now may be optimally 
leveraged to improve optical micromanipulation and the 
performance of light-driven machines. 
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